Simultaneous separation for the Kowalevski and 
Goryachev-C haply gin gyrostats 



Vadim B. KuznetsovQ 



Dept of Applied Maths, University of Leeds, LEEDS LS2 9JT, United Kingdom 

E-mail: vadim@maths.leeds.ac.uk 



Abstract 

In the special case of zero square integral the Kowalevski gyrostat and Goryachev-Chaplygin 
gyrostat share a simple separation of variables originated from the 4x4 Lax matrix. 



^EPSRC Advanced Research Fellow 



1 



Introduction 



The theory of separation of variables (SoV) is the theory of special canonical transformations 
and the theory of quantum separation of variables is the theory of corresponding integral 
transforms. The former was already understood in this sense starting from the development 
of the Hamilton- Jacobi approach for solving Liouville integrable systems. The real power of 
the latter, which by now is undoubted, requires further investigation and demonstration. 

The speciality of a separating transformation stems from its definition as a transform 
resulting in new variables being separated or equations being decoupled from one another. 
Below we give a (working) definition of SoV in the context of finite-dimensional integrable 
Hamiltonian dynamics. 

By separation of variables for an n-degrees-of-freedom integrable system having n inde- 
pendent Poisson commuting integrals {i/j (q, p)}"^]^ we mean a canonical transformation 
from the (old) Darboux variables qj,Pj, j = to new Darboux variables Uj,Vj, 

j = 1, . . . ,n, which satisfy the following separation equations: 



n 

E 



aij{ui,Vi) Hj{q,p) = bi{ui,Vi), i = l,...,n. (1) 



In other words, being expressed in terms of the new variables, the integrals of motion Hj 
acquire the following 'separated form': 

H = A-iB, (A)^^. = a,„ (B), = 6„ (H), = (2) 

The conditions that the functions aij and bi in (|I|) depend on the new (separation) variables 
with the index i only, is crucial. It indeed means that the n equations in (|^) are really 
separated from one another. Therefore, they are n equations, each of one degree of freedom, 
sharing only the common values of the Hamiltonians Hj 

The above definition includes, as even more special canonical transform, the (classical) 
coordinate separation of variables when the new coordinates, say u, are the functions of the 



old coordinates (q) only, and they do not depend on the momenta (p). See |jT8|, |T0| for 

some history of this sub-class of transformations and for many examples of such situation. 

General separating canonical transforms, however, have new (separation) variables which 
are non-trivial functions of all 2n initial canonical variables, 

Uj = Uj{ci,p), Vj=Vj{ci,p), j = l,...,n, (3) 

and which satisfy the standard Poisson brackets 

{uj,Uk} = {vj,Vk} = {uj,Vk} = 0, j ^ k, {vj,Uj} = l, j = l,...,n. (4) 

Examples of such separating canonical transforms are usually much more sophisticated than 
those of the coordinate ones. As far as I know, the first explicit example was given by van 
Moerbeke in 1976 in concerning the separation of variables for the periodic Toda lattice 



see also [|TT1)- In 1980 Kozlov rewrote the classical results of Goryachev and Chaplygin 



on integration in quadratures of the Goryachev- Chaplygin top as a simple canonical trans- 
form. The method of SoV, viewed as a method of special canonical transformations, was 

^as one can see, this definition is similar to the one defining Stackel systems 
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further developed by Komarov in a series of works on tops, including quantum separation of 
variables, see Many further examples have been produced since 1982, with the 

theory benefiting mostly from the developments of the algebraic geometric and r-matrix un- 
derstanding of the method of separation of variables. This led to a rather satisfying picture 
of the present state-of-art of non-coordinate separation of variables. See |Q and [jl6|, ^ for 
more details. 

Notice here that in the realm of the algebraic geometry usually associated with many 
Liouville integrable systems, the separation equations (|l]) appear as equations for a set of 
separating Darboux coordinates on the spectral curve of a Lax matrix L{u), namely: 

(1) ^ det{L{ui) - Vi) = 0, i = l,...,n. (5) 

SoV is not unique, so that in the situation when there exist several Lax matrices for the 
same integrable system, one should expect several separations. But even for the same Lax 
matrix there are many different separations given by different sections 0, cf. |]36| , p9| , ^ |. 



In the present paper we find a new canonical transform which simultaneously separates 



two classical tops: the Kowalevski top and the Goryachev-Chaplygin top ||T^, The 
latter top is integrable only for the zero value of the square integral, £ = (which is one of 
the Casimirs), so that although the canonical transform will be defined for arbitrary value 
of £, it will separate the Kowalevski top (and the other one) only for i = 0. SoV for the 
Kowalevski top for i 0, which would generalize the found transform, remains an unsolved 
problem. 



Sretenskii |p5[ discovered an integrable extension of the Goryachev-Chaplygin top adding 



the gyrostatic term to the Hamiltonian. Komarov |]2J] and, independently, Yehia fS^j found 



the gyrostat extension of the Kowalevski top in 1987. The quantum Goryachev-Chaplygin 



gyrostat was treated in . For quantum Kowalevski gyrostat see |2J] . Here we will always 
include such gyrostatic terms (cf. the parameter c below). 

In Section 1 we give definitions of the integrable systems concerned and their Lax matri- 
ces. In Section 2 we recall the separation method with main results presented in Section 3. In 
Section 4 we re-write the found separating transformation through the generating function. 
Solution of the inverse problem is given in Section 5. Finally, some concluding remarks can 
be found in the last Section. 



1 Algebra e(3), tops and Lax matrices 

The Poisson brackets for the e(3) generators Jk, Xk, k = 1,2,3, are defined in the standard 
way: 

{Jk, Ji} {xfc,Xi} = 0, (6) 

where Skim is the completely anti-symmetric tensor, £123 = 1. 
The Casimirs of the bracket (^) have the form 

Ci = xiJi + X2J2 + X'iJ'i = i, C2 = x\ + x\ + xl = 1. (7) 

The Kowalevski gyrostat has the Hamiltonian H, 

H = + {J3 + cf -2bxi, (8) 
^those usually differ by the number of non-moving poles of the Baker-Akhiezer function 



3 



and the second integral K, 



K={J3 + cfJ^ + 26(J3 + c)(JiX3 - -hxi) - b^l - 2hUu (9) 
which are Poisson commuting: 

{H,K] = Q. (10) 
Here we use the following notation for the square of the vector of angular momentum: 

j':=Jl + Jl + Jl (11) 

The two integrals of motion, H and K, define what is called Kowalevski gyrostat. It is a 
Liouville integrable system with two degrees of freedom. 

Another integrable system closely related to the Kowalevski gyrostat is the Goryachev- 
Chaplygin gyrostat. It is integrable only when £ = and is defined by two Poisson commuting 
integrals of motion: 



H = Jl + Jl + i2J,, + cf -Abxi, 
k = {J3 + c){Jl + Jl) + 2hx3Ji. 



(12) 
(13) 



The Kowalevski top without the gyrostatic term (c = 0) was integrated by Kowalevski 
in 1889 p9| . See also |^ where a 2 x 2 Lax matrix was constructed for the Kowalevski top 
which is related to the separation of variables hidden in Kowalevski's integration. There is a 
large body of literature dedicated to the Kowalevski top, including the study of its geometry. 
See |2g, 1^, ^, ig, 0, §, |TB|, §, ^ 0, ||, |§ to name just a few references. 

The Goryachev-Chaplygin top and gyrostat also received much attention in the literature 
starting from the discovery of this integrable top by Goryachev and Chaplygin in 1900. We 

I, 0, |,[ll, 0,13. 



refer the reader to the works 13, 0, 30, 21, 22 



In parallel with the notations Ji, J2 and Xi, X2 we will be using their equivalent complex 
versions J+, J_ and x+, X-. 



J± = Ji ± iJ2, x± = Xi ± 1x2, 

with the e(3) Poisson brackets (^ replaced, correspondingly, by 

{-^3, -^±} = Ti-^±, {J+, J-} = -2iJ3, {J3,a;±} = {a^3, = =Fia;±, 

{ J+, x_} = J_} = -2ix3, { J3, X3} = { J+, x+} = { J_, x_} = 0, 

{xk,xi] = Q, kj = ±,3, 

and the Casimirs (|^) by 

2Ci = x+ J_ + x_ J+ + 2x3 J3 = 2i, C2 = x+x. + X3 = 1. 



(14) 

(15) 

(16) 
(17) 



A 4 X 4 Lax matrix for the Kowalevski gyrostat was found in and used in ||^ to 
integrate the problem in terms of Prymian theta-functions. This integration is different 
from the one performed by Kowalevski. It is not obvious how to construct a separation of 
variables related to such Lax matrix. With the general case still being an interesting and 
challenging problem, we show here how to do it in the special case, when i = 0. That is, 
for this special case we construct a new separation of variables for the Kowalevski top (and 
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simultaneously for the Goryachev-Chaplygin top) with separation variables belonging to the 
spectral curve of the 4x4 Lax matrix from |5T[| . 
The required Lax matrix is as follows: 



L{u) = -i 



/ - 

u 
bx+ 

V? 

J+ 
u 
bx3 



bx3 

7(2 



J_ 

u 
bx-s 

V? 

2J3+C 



6x3 \ 



\ 0x3 J__ 



-2 - 



bx- 



2 + ^ 



(19) 



The spectral curve F: det(L(M) — t^) = 0, of the Lax matrix ([191) has the form 



r : 



+ — -— -4 



or, 



„/ H 4c^ H^-AK 262(2£2_m 
2t;2 2 - — + — - — + + i + ^ = 0. 



(20) 



(21) 



As I noticed in |^ (cf. also ||3]) this Lax matrix contains the 3x3 Lax matrix L{u) for 
the Goryachev-Chaplygin top as its (1, l)-minor: 



^ _£ 6x3 



L{u) = -i 



bx3 2J3+C 
u 

\ r,, n,Z 



u 

2 + ^ 

2J3+C 



:22) 



The spectral curve V: det(L(M) — w) = 0, of the Lax matrix (p^ has the form 
'2 / H b'^\ 4ic . .2k -cH .,^c^2x^l 



cv 



i 4_ + U + — + i- 



0. 



u 



u 



(23) 



The Lax matrix (|22D (and the curve (^3])), for i = 0, was used in to construct explicit 
theta-function formulas solving the dynamics of the Goryachev-Chaplygin top. 

In the present paper I take one step further in studying a close connection between two 
problems and show that there exist separation variables Uj,Vj, j = 1,2, which are canonical 
and which, for i = 0, belong to both curves, F and f , simultaneously: 



det{L{uj) — Vj) = 0, 
det{L{uj) — Vj) = 0, 



J = 1,2, 



{uj, Uk} = {vj, Vk} = {uj, Vk} = 0, j k, {vj, Uj} 



J 



1,2. 



(24) 



(25) 



Therefore, I construct a new separation of variables which is characterized by the property 
of being a simultaneous separation for both tops. 



2 The method of SoV 

Separation variables had been generally used to construct closed expressions for the action 
variables (in terms of abelian integrals) or to get a separated representation for the action 
function. Therefore, the SoV method had served for a long time an important but technical 
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role in solving Liouville integrable systems of classical mechanics. A new and much more 
exciting application of the method came with the development of quantum integrable sys- 
tems. Because of the fact that quantization of the action variables seemed to be a rather 
formidable task, quantum separation of variables became an inevitable refuge. In fact, it 
has been successfully performed for many families of integrable systems (see, for instance, 
survey |5^). 

Starting from about 1982 the method of separation of variables gets connected with the 
i?-matrix formalism of the quantum inverse scattering method, developed during that time 
by the Leningrad School. It was noticed by Komarov (see and |Q) that for the 2x2 
L-operators (Lax matrices) the separation variables ought to be the zeros of the off-diagonal 
element of the L-operator. This observation was fully exploited by Sklyanin in [^, |53 



who developed a beautiful (pure algebraic) setting for the method within the framework 
of the i?-matrix technique. Since then this approach took off and led to separations for 
many families of integrable systems. The method was further generalized to include higher 



rank L-operators and non-standard normalizations, see the 1995 review and the later 
developments in ||, g^, 0, |D|, ig. 

An alternative, algebraic geometric approach, which dates back to Adler and van Moer- 
beke |^ ^ and Mumford and includes many researchers, have been developed starting 
from about the same time (see, for instance, |15|, |l], ||, |5^, |1^, |T^). It is based on thorough 
studies of the geometry of lower genus (algebraic completely) integrable systems. It has also 
led to many important new separations for complicated systems and tops. 

Below we recall the most important formulas of the SoV method, adopting the algebraic 

mm)- 



description and following mainly the work 



(see also 



A Baker- Akhiezer function / is the eigenvector of the Lax matrix for an integrable system 



Hu) f = vf, 



(26) 



considered as a function on the spectral curve F : det{L{u) — v) = 0. The inverse scattering 
method pins down the separation variables as poles of this function (see the first exam- 
ple of this general fact in and |11[). There is, however, a large freedom of similarity 
transformations of the Lax matrix, 



Liu) 



VL{u)V- 



(27) 



which do not change the spectrum of L{u) but change the divisor of poles of /. This 
freedom can be characterized, and therefore fixed, by introducing a normalization of the 
Baker-Akhiezer function, 



N 



^ - f = ^OLi fi = l 



f = (/li • • • > / 



N 



[2S) 



which is given by a normalization (row-) vector a = (ai, . . . , Oat). In other words, one 
considers a section of the line-bundle (cf. \i7\ , 0). A proper (separating) normaliza- 
tion/section should give a divisor V = Yl'j=i{'^jiVj) of moving poles of /, consisting of n 
(independent) points on the curve F whose coordinates are canonical variables. The canon- 
icity of the separation variables is usually checked by a calculation involving a r-matrix, 
but in lower genus situations it can be proved in a direct calculation. For a non-separating 
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normalization/section the divisor V will usually have more points than needed, which will 
not give canonical variables. 

Many families of Lax matrices have a simple separating normalization when one of the 
components of the Baker-Akhiezer vector is put 1 (and, hence, one looks at the common 
poles of the other components), which corresponds to the following vector 5 0: 

ao = (1,0,...,0,0). (29) 

We will call such normalization the standard normalization. There are examples of non- 



standard (dynamical) separating normalizations for the systems with elliptic r-matrix 



T7|| , Calogero-Moser systems and the D-type Toda lattice [Q. See also where 
non-standard dynamical separating normalizations were used to construct Backlund trans- 
formations. Generally, Lax matrices with extra symmetries require non-standard separating 
normalizations. 

Now, assuming that we know a separating normalization a for an integrable system 
with the Lax matrix L(m), let us derive the equations for the separation variables {uj,Vj), 
j = l,...,n. 



From the linear problem (^) and normalization 
0, . . . , — 1, hence. 



/ 



a 

a ■ L{u) 

N-li 



\ 



\ a-L^-\u) J 



we derive that a ■ L f 
1 

V 

N-l 



k 



(30) 



Another useful representation for the eigenvector /, which can be verified directly, is as 
follows: 

{L{u) - 
fj — 



(a 



L(u)-v) 



yk = l,...,N, 



(31) 



where the wedge denotes the adjoint matrix. It follows from it that the poles of / are the 
common zeros of the vector equation: 



a ■ {L{u) — v)' 



0. 



(32) 



Eliminating v from these equations, one can get a single equation for the w-components of 
the separating variables as zeros of the following determinant: 



B{u) = det 



d ■ L{u) 



0. 



(33) 



Notice that this is exactly the denominator in the representation (^). 

Also, from the equations (|3^ ) we can obtain explicit formulas for the u-components of 
the separation variables in the form 

v = A{u), (34) 
with A[u) being rational functions of the entries of Liu). Let us derive those formulas, 
^or to any other similar vector with the 1 elsewhere 
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Define the matrices L^p\ p = 1, . . . , N, with the following entries: 



N N 
h=l ip_i=l 



^5^p— 1 



1'p—l ;^p— 1 



, p = 2,3,...,iV, 



(35) 



and put L^^^ = L. These matrices satisfy the recursion relation of the form 

L^p) = L (tr L(P-i)) - (p - 1) L^-^^L . 
Introduce the matrix B{u) by the formula 



(36) 



B{u) := 



I a- L('\u) L~\u) \ 
a-L^^\u) L~\u) 



\j^ a-LWi^u) L-\u)j 



With the help of this matrix we can represent the system of equations 
linear homogeneous equations for the vector of powers of (— f): 



a.[L{u)-v)^^{{-vr-\{-vr-\ 

from which we derive that 

{B''{u))k^ 



1) ■ B{u) = 0, 



VA;. 



(37) 

as a system of 
(38) 

(39) 



The formulas (|39D give plenty of representations for the rational functions A{u) in (0), all 
of them being compatible on the separation variables since, because of the equality 



Biu) 



det{B{u)) 



(40) 



the matrix B^{uj) has rank 1. For more details see 



3 A new separation of variables 

Consider the standard normalization vector 

ao= (1,0,0,0) 

for the 4x4 Lax matrix ([19|) . Then, for i = 0, the defining equations. 



im - v] 



Ik 



0, 



A; = l,2,3,4, 



give the following polynomial B{u) for the separation variables Ui and U2'- 

B{u) =u^ + B2U^ + Bo = {u^ - ul){u^ - ul), 



(41) 



(42) 



(43) 
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2b{Js + c)(x3J- - X-J3) + b^{xt + 2x1) 2b^X3X-{cJ- + bx. 



_ 2 ((J3 + C)J-+&X3)' 

J2(j! + 2te_) • ^^^^ 
Note also the following useful formula for the variable B2: 



^2 = -12^. (:^-^-^3 - (J' + bx^)x,) , (46) 



which is a linear expression in terms of J3 and 0:3. 

The corresponding rational function A{u) can be chosen as follows: 

A{u) = ^ + ^ , (47) 

. , iX3(J!+2bX-) .^ (J3 + C)J_+&X3 

A_i = icH , A_3 = 16 . (48) 

(V (V 

Now, it is easy to check the following Poisson brackets between the two functions: 

{A{u),A{v)} = {B{u),B{v)} = 

^u,veC. (49) 

{A{u),B{v)} = ^j,^ {u'B{v) - v'B{u)) , 

Using these formulas one derives that the variables Ui, U2, defined by (^3|)-(^6|), and their 
conjugated counterparts, 

Vj = A{uj), 3 = 1,2, (50) 
are indeed canonical (Darboux) variables: 

{uj,Uk} = {vj,Vk} = {uj,Vk} = 0, j 7^ k, {vj,Uj} = l, i = l,2. (51) 

These variables by construction satisfy (for £ = 0) to ( ^2]) and, therefore, to (|2^) . 



Notice here that the definitions of the separation variables do not depend on the value 
of £ and also that the brackets (pl|) are true for any i. 



4 Generating function 

Let us fix a special representation of the underlying e(3) algebra which will allow us to write 
down the found canonical transformation explicitly through the generating function. It is a 
kind of holomorphic representation in terms of the Darboux variables qj, pj, j = 1,2:0 

J_ = gi, x_ = q2, (52) 

J3 = i{qipi + q2P2) +£, X3 = igaPi + 1, (53) 
J+ = qipi + 2q2PiP2~2iipi-2ip2, x+ = q2pl - 2ipi. (54) 
"'cf. a holomorphic representation of sl(2) 
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The meaning of this reahzation is in the fact that the variables J_ and X- do not depend 
on the momenta p and the variables J3 and .T3 arc linear in momenta. This together with 
linearity of the variables U1U2 and u\ + u\ in terms of J3 and make it possible to integrate 
the equations explicitly in terms of elementary functions. 

Finally, the canonical transformation defined in the previous Section is given by the 
following generating function F(u|q): 

F(u|q) = ^ + \og{ql + ^hq^) + mogiq^) + 



^2 ^ 2uiU2Jqi + 2bq2 



+ iy^FT^f2^-^). (55) 
2 ^ \ b q2UiU2 J 

Therefore, the equations of the change of variables (q, p) <-> (u, v) are written in the form 

5 Inverse problem 

The inverse problem, i.e. finding expressions for the initial e(3) variables in terms of sepa- 
ration variables, also has an explicit solution which is given below: 

2 1,2 '^l('^2'^2 ~b ^ — li2('^l'^l ~ b ) ( 2/ 4 2 7,2\ 2/ 4 2 7,2\ , ,1 2 2/ 2 2^ 

?1 = 7-^-6? X2 («1 %^2 -b)- M2 - & + 4Mi«2 - 

'm\'u!2\U\V\ — U2V2) ^ 



(57) 

6(^.2-^2) (^2(^4^2 _ ^2) _ ^2(^4^2 _ ^2)) 
J_ — qi, — n 4 A/ \2 ' '^^^ 



2uiul{uivi - U2V2) {ul{u\vl - 62) - ul{u\vl - 62)) 
X \b'^{u\ - ulf + 2cu\ul{uiVi - U2V2) {u\{yl - 2) - ■U2('i^2 ~ 2)) 

+ 2dP'u\u\{uiV\ — U2V2){u\ — ul) + 2ib^ulul{ulvl — ulv2){ul — u^) 
+ iUiU^iuiVi - U2V2) {viu\{vl - 4) - V2u\{vl - 4) + uiU2ViV2{uiVi - U2V2)) ) , (59) 

_ 2igiu2^2 (^ul{uiVi - ic) - ul{u2V2 - ic) 

6(?/2(^xt?;22 - 62)- 7/2(^,4^2 _ ^2)) ' l^^j 

J 4:qiu\ul{uiVi - U2V2) ^ 

^ ~ 62 {uliuivi - 62) - uliufv! -h^)f ^ 

X ^ u\u2ViV2{uiVi — U2V2) + b'^{u\vi — U2V2) — idP'{u\ — ul) 

— iculuKulvf — u\v2) — (?u\u\{uiVi — U2V2) ) , (61) 

_ 1 

~ 6 {ul{u\vl - 62) - ul{u\vl - ^ 
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X -Aicy^ujuKul - ul)(ulvi - ulv2) + 2b'^ujul(uj - ul)(ufvj - U2vl) 
-2b^c\juliuj - ulf - 2c^u\u\{u\ - ul) [u\{v\ - 4) - ul{vl - 4)) 
- ^\cu\u\{u\v^ - ulvi) (u\{y\ - 4) - u\{yl - 4)) 

+ 2u\u\ ({u\v\ - ulvl) {u\vl - uy^) - A{u\v^ - ulv^f) ) . (62) 

6 Concluding remarks 

A separation of variables for the Kowalevski gyrostat, with £ = 0, is found for the first time. 
It appeared to separate the Goryachev-Chaplygin gyrostat as well, thereby giving another 
separation for this problem. 

If we put c = 0, i.e. switch off the gyrostatic term, we can compare our results with 
the original Kowalevski's and Goryachev-Chaplygin 's separations for the respective tops (in 
the £ = case). It is easy to see that the new separation is as complicated as the (simple) 
Goryachev-Chaplygin separation and it is much simpler than the Kowalevski separation. 

Hence, the new separation stands a good chance to be quantized. 
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